What Is Claimed Is: 



1 LA method for using a computer system to solve an unconstrained 

2 interval global optimization problem specified by a function / wherein/is a scalar 

3 function of a vector x = {x h x 2 > x 3) ... x n \ the method comprising: 

4 receiving a representation of the function / at the computer system; 

5 storing the representation in a memory within the computer system; and 

6 performing an interval global optimization process to compute guaranteed 

7 bounds on a globally minimum value of the function^*) over a subbox X; 

8 wherein performing the interval global optimization process involves, 

9 applying term consistency to a set of relations associated 

1 0 with the function / over the subbox X, and excluding any portion of 

1 1 the subbox X that violates any of these relations, 

12 applying box consistency to the set of relations associated 

1 3 with the function/ over the subbox X, and excluding any portion of 

14 the subbox X that violates any of these relations, and 

15 performing an interval Newton step on the subbox X to 

1 6 produce a resulting subbox Y, wherein the point of expansion of 

1 7 the interval Newton step is a point x within the subbox X, and 

1 8 wherein performing the interval Newton step involves evaluating 

19 the gradient g(x) of the function/x) using interval arithmetic. 

1 2. The method of claim 1 , wherein applying term consistency 

2 involves: 

3 symbolically manipulating an equation within the computer system to 

4 solve for a term g(xj), thereby producing a modified equation g(x 3 ) = h(x\ wherein 

5 the term g(Xj) can be analytically inverted to produce an inverse function g^iy); 
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substituting the subbox X into the modified equation to produce the 
equation g(X)) = h(X); 

solving for X) = g\h(X)); and 

intersecting X) with the interval Xj to produce a new subbox X + ; 

wherein the new subbox X + contains all solutions of the equation within 
the subbox X, and wherein the size of the new subbox X + is less than or equal to 
the size of the subbox X. 

3. The method of claim 1, wherein performing the interval global 
optimization process involves: 

keeping track of a smallest upper bound f_bar of the function/x); 
removing from consideration any subbox X for which/X) >f_bar; and 
wherein applying term consistency to thenar relation involves applying 
term consistency to the inequality X x ) * f_bar over the subbox X. 

4. The method of claim 3, wherein applying box consistency to the 
set of relations involves applying box consistency to the inequality /x) <f_bar 
over the subbox X. 

5. The method of claim 1, wherein performing the interval global 
optimization process involves: 

determining the gradient g(x) of the function fix), wherein g(x) includes 
components g,( x ) (* = 1, •■■>«); 

removing from consideration any subbox for which any element of g(x) is 
bounded away from zero, thereby indicating that the subbox does not include a 
stationary point of fix); and 
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8 wherein applying term consistency to the set of relations involves applying 

9 term consistency to each component g,(x)=0 (i=l, ...,«) of g(x)=0 over the subbox 
10 X. 

1 6. The method of claim 5, wherein applying box consistency to the 

2 set of relations involves applying box consistency to each component 

3 g,(x)=0 (i=l, ...,«) of g(x)=0 over the subbox X. 

1 7. The method of claim 1 , wherein performing the interval global 

2 optimization process involves: 

3 determining diagonal elements H n (x) (7=1, ...,«) of the Hessian of the 

4 function^*); 

5 removing from consideration any subbox for which a diagonal element of 

6 the Hessian is always negative, which indicates that the function/is not convex 

7 and consequently does not contain a global minimum within the subbox; 

8 wherein applying term consistency to the set of relations involves applying 

9 term consistency to each inequality H n {x) > 0 (z=l, ...,«) over the subbox X. 

1 8. The method of claim 7, wherein applying box consistency to the 

2 set of relations involves applying box consistency to each inequality 

3 H u (x) > 0 (z=l, ...,«) over the subbox X. 

1 9. The method of claim 1, 

2 wherein performing the interval Newton step involves, 

3 computing the Jacobian J(x 5 X) of the gradient g evaluated 

4 as a function of a point x over the subbox X, 
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computing an approximate inverse B of the center of 
J(x,X), and 

using the approximate inverse B to analytically determine 
the system Bg(x) 5 wherein g(x) is the gradient of the function/x), 
and wherein g(x) includes components g l (x) (i=l, ...,«); and 
wherein applying term consistency to the set of relations involves applying 
term consistency to each component (Bg (x)), = 0 (f=l, ...,«) for each variable 
x t (7=1, ...,«) over the subbox X. 

1 0. The method of claim 9, wherein applying box consistency to the 
set of relations involves applying box consistency to each component (Bg(x)), = 0 
(i=l, „. } ri) for each variable x t (i=l, ...,«) over the subbox X. 

1 1 . The method of claim 1 , further comprising terminating attempts to 
further reduce the subbox X when: 

the width of X is less than a first threshold value; and 

the magnitude of/(X) is less than a second threshold value. 

1 2. The method of claim 1 1 , wherein performing the interval Newton 
step involves: 

computing J(x,X), wherein J(x,X) is the Jacobian of the function f 
evaluated as a function of x over the subbox X; and 

determining if J(x 5 X) is regular as a byproduct of solving for the subbox Y 
that contains values of y that satisfy M(x,X)(y-x) = r(x), where 
M(x,X) = BJ(x,X), r(x) = -Bf(x), and B is an approximate inverse of the center of 
J(x,X). 
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13. A computer-readable storage medium storing instructions that 
when executed by a computer cause the computer to perform a method for using a 
computer system to solve an unconstrained interval global optimization problem 
specified by a function / wherein /is a scalar function of a vector 

x = (xj, X2, X3, ... x n ), the method comprising: 

receiving a representation of the function / at the computer system; 
storing the representation in a memory within the computer system; and 
performing an interval global optimization process to compute guaranteed 
bounds on a globally minimum value of the function^*) over a subbox X; 
wherein performing the interval global optimization process involves, 
applying term consistency to a set of relations associated 
with the function / over the subbox X, and excluding any portion of 
the subbox X that violates any of these relations, 

applying box consistency to the set of relations associated 
with the function / over the subbox X, and excluding any portion of 
the subbox X that violates any of these relations, and 

performing an interval Newton step on the subbox X to 
produce a resulting subbox Y, wherein the point of expansion of 
the interval Newton step is a point x within the subbox X, and 
wherein performing the interval Newton step involves evaluating 
the gradient g(x) of the function/x) using interval arithmetic. 

14. The computer-readable storage medium of claim 13, wherein 
applying term consistency involves: 

symbolically manipulating an equation within the computer system to 
solve for a term g(x y ), thereby producing a modified equation g(Xj) = h(x), wherein 
the term g(xj) can be analytically inverted to produce an inverse function g' J (y); 
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6 substituting the subbox X into the modified equation to produce the 

7 equation g(X , J ) = h(X); 

8 solving for X) = g\h(X)); and 

9 intersecting X) with the interval Xj to produce a new subbox X + ; 

10 wherein the new subbox X + contains all solutions of the equation within 

1 1 the subbox X, and wherein the size of the new subbox X + is less than or equal to 

12 the size of the subbox X. 

1 15. The computer-readable storage medium of claim 1 3 , wherein 

2 performing the interval global optimization process involves: 

3 keeping track of a smallest upper bound f_bar of the function^x); 

4 removing from consideration any subbox X for which/X) >f_bar; and 

5 wherein applying term consistency to the f_bar relation involves applying 

6 term consistency to the inequality fix) < f_bar over the subbox X. 

1 1 6. The computer-readable storage medium of claim 15, wherein 

2 applying box consistency to the set of relations involves applying box consistency 

3 to the inequality fix) < f_bar over the subbox X. 

1 17. The computer-readable storage medium of claim 1 3, wherein 

2 performing the interval global optimization process involves: 

3 determining the gradient g(x) of the function fix), wherein g(x) includes 

4 components g,(x) (f=l, ...,«); 

5 removing from consideration any subbox for which any element of g(x) is 

6 bounded away from zero, thereby indicating that the subbox does not include a 

7 stationary point of/(x); and 



31 

Attorney Docket No. SUN-P7077-SPL Inventors: Walster et al. 

ARPH \SUN MICROSYSTEMS\SUN-P7077-SPL\SUYN-P7077 APPLICATIONS DOC 



wherein applying term consistency to the set of relations involves applying 
term consistency to each component g,(x)=0 (i=l, of g(x)=0 over the subbox 
X. 

1 8. The computer-readable storage medium of claim 17, wherein 
applying box consistency to the set of relations involves applying box consistency 
to each component g,( x ) =0 (i=l, ...,«) of g(x)=0 over the subbox X. 

19. The computer-readable storage medium of claim 13, wherein 
performing the interval global optimization process involves: 

determining diagonal elements H n (x) (7=1, ...,«) of the Hessian of the 
function^x); 

removing from consideration any subbox for which a diagonal element of 
the Hessian is always negative, which indicates that the function /is not convex 
and consequently does not contain a global minimum within the subbox; 

wherein applying term consistency to the set of relations involves applying 
term consistency to each inequality H n (x) > 0 0=1, ...,«) over the subbox X. 

20. The computer-readable storage medium of claim 19, wherein 
applying box consistency to the set of relations involves applying box consistency 
to each inequality H n (x) > 0 0=1, ...,n) over the subbox X. 

2 1 . The computer-readable storage medium of claim 1 3, 
wherein performing the interval Newton step involves, 

computing the Jacobian J(x,X) of the gradient g evaluated 
as a function of a point x over the subbox X, 
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5 computing an approximate inverse B of the center of 

6 J(x,X), and 

7 using the approximate inverse B to analytically determine 

8 the system Bg(x), wherein g(x) is the gradient of the function/x), 

9 and wherein g(x) includes components g,(x) (i=l, ...,«); and 

1 0 wherein applying term consistency to the set of relations involves applying 

1 1 term consistency to each component (Bg (x)), = 0 (7=1, ...,ri) for each variable 

12 x t (7=1, ...,ri) over the subbox X. 

1 22. The computer-readable storage medium of claim 21, wherein 

2 applying box consistency to the set of relations involves applying box consistency 

3 to each component (Bg (x)), = 0 (7=1, ...,ri) for each variable x t (7=1, ... f ri) over the 

4 subbox X. 

1 23. The computer-readable storage medium of claim 13, wherein the 

2 method further comprises terminating attempts to further reduce the subbox X 

3 when: 

4 the width of X is less than a first threshold value; and 

5 the magnitude of/[X) is less than a second threshold value. 

1 24. The computer-readable storage medium of claim 13, wherein 

2 performing the interval Newton step involves: 

3 computing J(x,X), wherein J(x,X) is the Jacobian of the function f 

4 evaluated as a function of x over the subbox X; and 

5 determining if J(x,X) is regular as a byproduct of solving for the subbox Y 

6 that contains values of y that satisfy M(x,X)(y-x) = r(x), where 
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1 M(x 5 X) == BJ(x,X), r(x) = -Bf(x), and B is an approximate inverse of the center of 

2 J(x,X). 

1 25. An apparatus that solves an unconstrained interval global 

2 optimization problem specified by a function / wherein /is a scalar function of a 

3 vector x = (xj, x 2) X3, ... x n ), the apparatus comprising: 

4 a receiving mechanism that is configured to receive a representation of the 

5 function f\ 

6 a memory for storing the representation; and 

7 an interval global optimization mechanism that is configured to perform 

8 an interval global optimization process to compute guaranteed bounds on a 

9 globally minimum value of the function/x) over a subbox X; 

10 a term consistency mechanism within the interval global optimization 

1 1 mechanism that is configured to apply term consistency to a set of relations 

12 associated with the function / over the subbox X ? and to exclude any portion of the 

1 3 subbox X that violates any of these relations; 

14 a box consistency mechanism within the interval global optimization 

1 5 mechanism that is configured to apply box consistency to the set of relations 

16 associated with the function / over the subbox X 5 and to exclude any portion of the 

1 7 subbox X that violates any of these relations; and 

1 8 an interval Newton mechanism within the interval global optimization 

1 9 mechanism that is configured to perform an interval Newton step on the subbox X 

20 to produce a resulting subbox Y, wherein the point of expansion of the interval 

21 Newton step is a point x within the subbox X, and wherein performing the interval 

22 Newton step involves evaluating the gradient g(x) of the function/*) using 

23 interval arithmetic. 
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26. The apparatus of claim 25, wherein the term consistency 
mechanism is configured to: 

symbolically manipulate an equation to solve for a term g(Xj\ thereby 
producing a modified equation g(xj) = /z(x), wherein the term g(xj) can be 
analytically inverted to produce an inverse function g' l (y); 

substitute the subbox X into the modified equation to produce the equation 
g(X'j) = KX); 

solve for X) = g\h(X))- and to 

intersect X) with the interval Xj to produce a new subbox X + ; 

wherein the new subbox X + contains all solutions of the equation within 
the subbox X, and wherein the size of the new subbox X + is less than or equal to 
the size of the subbox X. 

27. The apparatus of claim 25, 

wherein the interval global optimization mechanism is configured to, 

keep track of a smallest upper bound f_bar of the function 
fix), and to 

remove from consideration any subbox X for which 
XX) >f_bar\ and 
wherein the term consistency mechanism is configured to apply term 
consistency to the inequality fix) < f_bar over the subbox X. 

28. The apparatus of claim 27, wherein the box consistency 
mechanism is configured to apply box consistency to the inequality fix) < f_bar 
over the subbox X. 

29. The apparatus of claim 25, 
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2 wherein the interval global optimization mechanism is configured to, 

3 determine the gradient g(x) of the function f[x\ wherein 

4 g(x) includes components g t (x) (z=l, and to 

5 remove from consideration any subbox for which any 

6 element of g(x) is bounded away from zero, thereby indicating that 

7 the subbox does not include a stationary point of fix); and 

8 wherein the term consistency mechanism is configured to apply term 

9 consistency to each component g,(x)=0 (zW, ...,ri) of g(x)=0 over the subbox X. 

1 30. The apparatus of claim 29, wherein the box consistency 

2 mechanism is configured to apply box consistency to each component 
D 3 &(x)=0 (f =7, ... } n) of g(x)=0 over the subbox X. 

m 

; ~ 1 31. The apparatus of claim 25 , 

4? 2 wherein the interval global optimization mechanism is configured to, 

3 determine diagonal elements H u (x) (/ = 1, ...,») of the 
l1 4 Hessian of the function /{x) 5 and to 

j2{ 5 remove from consideration any subbox for which a 

Q 6 diagonal element of the Hessian is always negative, which 

7 indicates that the function /is not convex and consequently does 

8 not contain a global minimum within the subbox; 

9 wherein the term consistency mechanism is configured to apply term 
10 consistency to each inequality H n (x) > 0 (/=1, ... } ri) over the subbox X. 

1 32. The apparatus of claim 3 1 , wherein the box consistency 

2 mechanism is configured to apply box consistency to each inequality 

3 H u (x) > 0 (z=l, ...,«) over the subbox X. 
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33. The apparatus of claim 25, 

wherein the interval Newton mechanism is configured to, 

compute the Jacobian J(x,X) of the gradient g evaluated as 
a function of a point x over the subbox X, 

compute an approximate inverse B of the center of J(x,X), 

and to 

use the approximate inverse B to analytically determine the 
system Bg(x), wherein g(x) is the gradient of the function/(x) 3 and 
wherein g(x) includes components g l (x) (i=l, ...,«); and 
wherein the term consistency mechanism is configured to apply term 
consistency to each component (Bg(x)), = 0 (7=1, ...,ri) for each variable 
(7=1, ...,ri) over the subbox X. 

34. The apparatus of claim 33, wherein the box consistency 
mechanism is configured to apply box consistency to each component 
(Bg(x))i = 0 (z=l, ...,«) for each variable x t (7=1, ...,«) over the subbox X. 

35. The apparatus of claim 25, further comprising a termination 
mechanism that is configured to terminate attempts to further reduce the subbox X 
when: 

the width of X is less than a first threshold value; and 

the magnitude of J(X) is less than a second threshold value. 

36. The apparatus of claim 11, wherein the interval Newton 
mechanism is configured to:, 
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1 compute J(x,X) 5 wherein J(x,X) is the Jacobian of the function f evaluated 

2 as a function of x over the subbox X; and to 

3 determine if J(x,X) is regular as a byproduct of solving for the subbox Y 

4 that contains values of y that satisfy M(x,X)(y-x) = r(x), where 

5 M(x,X) = BJ(x ? X) ? r(x) = -Bf(x), and B is an approximate inverse of the center of 

6 J(x,X). 
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